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ABSTRACT 

We present a clustering analysis of Luminous Red Galaxies (LRGs) in Stripe 82 from the 
Sloan Digital Sky Survey (SDSS). We study the angular two-point autocorrelation function, 
w{6), of a selected sample of over 130000 LRG candidates via colour-cut selections in izK 
with the K band coverage coming from UKIRT Infrared Deep Sky Survey (UKIDSS) LAS. 
We have used the cross-correlation technique of Newman (2008) to establish the redshift 
distribution of the LRGs. Cross-correlating them with SDSS quasi-stellar objects (QSOs), 
MegaZ-LRGs and DEEP2 galaxies, implies an average redshift of the LRGs to be z a; 1 
with space density, jtg ^ 3.20 + 0.16 x lO"'' h^'Mpc"^. For 6 < 10' (corresponding to a; 10 
h^'Mpc), the LRG w{6) significantly deviates from a conventional single power-law as noted 
by previous clustering studies of highly biased and luminous galaxies. A double power-law 
with a break at r/, a 2.4h"'Mpc fits the data better, with best-fit scale length, ro,i = 7.63 ± 0.27 
h^'Mpc and slope yi - 2.01 ± 0.02 at small scales and ro.2 - 9.92 + 0.40 h^'Mpc and 
yi-l .64 + 0.04 at large scales. Due to the flat slope at large scales, we find that a standard A 
cold dark matter (ACDM) linear model is accepted only at 2 - 3cr, with the best-fit bias factor, 
b - 2.1 A + 0.07. We also fitted the halo occupation distribution (HOD) models to compare our 
measurements with the predictions of the dark matter clustering. The effective halo mass of 
Stripe 82 LRGs is estimated as M^s - 3.3 ± 0.6 x 10'^ h 'Mq. But at large scales, the current 

HOD models did not help explain the power excess in th e clustering signa l. 

We then compare the w{ff) results to the results of ISawangwit et afl (1201 ll) from three 
samples of photometrically selected LRGs at lower redshifts to measure clustering evolution. 
We find that a long-lived model may be a poorer fit than at lower redshifts, although this as- 
sumes that the Stripe 82 LRGs are luminosity-matched to the AAQ. LRGs. We find stronger 
evidence for evolution in the form of the z ~ 1 LRG correlation function with the above flat 
2-halo slope maintaining to i > 50h 'Mpc. Applying the cross-correlation test of iRoss et al.l 
(l20Tlh . we find little evidence that the result is due to systematics. Otherwise it may repre- 
sent evidence for primordial non-Gaussianity in the density perturbations at early times, with 
= 90+30. 

Key words: galaxies: clustering - luminous red galaxies: general - cosmology: observations 
- large-scale structure of Universe. 



1 INTRODUCTION 

The statistical study of the clustering properties of massive galaxies 
provides important information about their formation and evolution 
which represent major questions for cosmology and astrophysics. 
The correlation function of galaxies remains a simple yet power- 
ful tool for im plementing such statistical clustering studies, (e.g. 
IPeeblesll 19801) . 

A lot of interest has been concentrated specifically on mea- 
suring the clustering correlation function of luminous red galax- 
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ies (LRGs) jEisenstein et al.l 2001 


) (see e.g IZehavi et al.l 2005bl; 


Blake, Collister. & Lahav 


l2008l: 


RossetalJ 20081: Wake et al.l 


20081; Sawanewit et alj 


2011). LRGs are credominantlv red 



massive e arly-type galaxies, intrinsically luminous (> 3L*) 
teisenstein et al. 2003; Loh & Strauss 2006; Wake et al. 200^ and 
thought to lie in the mo st massive dark matter halo es. They are also 
strongly biased objects jPadmanabhan et al.l2007l) and this coupled 
with their bright luminosity makes their clustering easy to detect 
out to high redshifts. For linear bias, the form of the LRG corre- 
lation function will trace that of the mass but even in this case the 
rate of correlation function evolution will depend on the bias model 
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(e.g. lFrvl 19961) , which in turn depends on the galaxy formation pro- 
cess. 

The passive evolution of the LR G LP and slow evoluti on of the 
LRG clustering jWake et al . 2008: Sawangwit et al1l201 ih seen in 
SDSS, 2SLAQ and AAQ. Surveys already presents a challenge for 
hierarchical models of galaxy formation as predicted for a cold dark 
matter (CDM) universe. Since the LRG clustering evolution with 
redshift has been controversial, a major goal is to use the angular 
congelation function to test if the slow clustering evolution trend 
continues out to z a 1. 

The uniformity of the LRG Spectral Energy Distributions 
(SEDs) with their 4000A Call H&K break, offer the ability 
to apply a colour-colour selection algorithm for our candidates. 
This technique has been successfully demonstrated primarily by 
lEisenstein et al.l in SDSS in the analysis of LRG cl ustering at 
low redshift and then in 2SLAQ jCannon et al]|2006h and AAQ. 
( IRos s et al. 200^ LRG surveys at higher redshifts. For our study, 
the available deep optical-IR iigrizJHK imaging data from the 
SDSS + UKIDSS LAS/DXS surveys in Stripe 82 will be used. This 
combination of NIR and deep optical imaging data, on a moderate 
sample size of area ~ 200 deg^, results in a sample of a 130000 
LRG candidates at redshift z ~ I- 

The main tool for our clustering analysis will be the two- 
point angular correlation function, w(d), which has been frequently 
used in the past, usually in cases where detailed redshift infor- 
mation was not known. Hence, selecting Stripe 82 LRGs based 
on colour-magnitude criteria, correspond to a rough photometric 
redshift (photo-z) estimation based on the 4000A break shifting 
through the passbands. We shall apply the cross- correlation tech- 
nique which was introduced by iNewmarJ ( l2008h to measure the 
redshift distribution, n(z), of our photometrically selected samples. 
One of the main advantages of w(6) is that it onl y needs the n (z) 
of the sample and then through Limber's formula dLimbedl 19531) it 
can be related to the spatial two-point correlation function, ^(r). 

In recent clustering studies, it was noted that the behaviour 
of f(r), which has previously been successfully described by a 
single power-law of the form ^(r) = (rlroy, significantly de- 
viates from such a power-law at ~ 1 h 'Mpc. The break in the 
power-law, can be interpreted in the framework of a halo model, 
as arising from the transition between small scales (1-halo term) to 
larger than a single halo scales (2-halo term). Currently, our the- 
oretical understanding of how galaxy clustering relates to the un- 
derlying dark matter is pro vided by the halo occupatio n distribu 



tion model (HOD, see, e .g Jing, Mo, & Boemer 1998 ; Ma & Fry 



2000| : jPeacock & Smitd |2000: Seliak .2000. : Scoccimarro et al.. 
200 ll : iBerlind & Weinberg! I2OO2I) via dark matter halo bias and 



halo mass function. Furthermore, the evolution of HOD can also 
give an insigh t into how certain galaxy populations evolve over 
cosmic ti me (White etal. 2007; S eo, Eisenstein. & Zehavill2008l : 
IWake et aU 2008 : Sawangwit et al..,2mlh . 

The outline of this paper is as follows. In Section 2, we briefly 
describe the SDSS and UKIDSS data used in this paper, while in 
Section 3 we describe the angular function correlation function es- 
timators and their statistical uncertainties. In Section 4, we estimate 
the redshift distribution through cross-correlations and then present 
the correlation results together with their power-law fits, ACDM 
model and a halo model in Section 5. Section 6 is devoted to in- 
teipretation of the clustering evolution. In section 7, we explore 
potential systematic errors that might affect the large scale clus- 
tering signal. We then argue that, if real, an observed large-scale 
clustering excess may be due to the scale-dependent bias caused by 
primordial non-Gaussianity and compare our results to other pre- 
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Figure 1. iz vs zK colour-colour plot. Priority A and B coiTcspond to the 
~ 700deg"^ and ~ 240deg"~ LRG samples, respectively. Objects with J - 
K < 1.3 which is typical for M stars are pfotted as bfack circfes where as 
those with J - K ^ 1.3 are pfotted in green. Evotutionary tracks for singfe 
burst (red fine) and r = IGyr (bfue fine) are overpfotted from z = to f .6 
with symbofs indicate z intervaf of 0.2. The evotutionary track of fate type 
gafaxies (magenta fine) is afso shown for comparison. 



vious works in Section 8. Finally, in Section 9 we summarize and 
conclude our findings. 

Throughout this paper, we use a flat A-dominated cosmology 
with U„ = 0.27, Ho = IOO/7 kms-'Mpc"', h=0.7, erg = 0.8 and 
magnitudes are given in the AB system unless otherwise stated. 



2 DATA 

2.1 LRG sample selection 

We perform a A^-band selection of high redshift LRGs in Stripe 
82 based on the co mbined optical and IR imaging data, ugrizJHK, 
from SDSS DR7 jAbazaiian et al]|2009l) and U KIDSS LAS sur- 
veys jLawrence et al.l 20071 : Warren et al.] l2007b . respectively. In 
previous studies, gri and riz colours have been used to select low 
to mediu m redshift LRGs, su ch as SDSS jEisenstein et al. I l200l[) . 
2SLAQ ( ICannonetaLll2006l) and AAQ. i IRoss etalj|2008h LRGs 
surveys up to z a 0.7. In this work we aim to study LRGs at z a 1, 
thus we use the izK colour magnitude limits for our selection in or- 
der to sample the 4000A Call H&K break of the LRGs' SED as it 
moves across the photometric filters (Fukugita et al. 1996: Smith et 
al. 2002) taking advantage of the NIR photometry coverage from 
UKIDSS LAS. Coupling the UKIDSS LAS to Kvega < 18 with the 
SDSS ugriz imaging to i^B < 22.5 in Stripe 82 produces an unri- 
valed combination of survey area and depth. Our selection criteria 
are : 



-2 



SDSS Best Imaging 

Z-K + 0.9{i -z)> 1.8, Pri A ~ lOOdeg 
Z-K + 0.9{i -z)>2.3, Pri B ~ 240deg-^ 
z-K- 0.9(i -z)> -0.2 
-0.5 <i-z< 1.7 
z-K< 4.0 
17.0 < a: < 18 

z < 22.0. 



(1) 
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The photometric selection of LRGs at z > 1 requires a combi- 
nation of optical and NIR photometry as the 4000A band straddles 
the z band. The selection of high-redshift LRGs is done on the ba- 
sis of SDSS iz photometric da ta and the LAS K band da ta (Fig.[T]l. 
LRG evolutionary models of ISruzual & Chario ^ ( l2003h are over- 
plotted for single burst and t = IGyr galaxy models indicating the 
izk plane area where we should apply our selections in order to 
study the high-z LRG candidates. 

Late-type star contamination is a major problem in selecting 
a photometric sample of z ~ I LRGs. Here the z - K colour also 
helps to distinguish the M stars colour locus from those of galaxies. 
From Fig[T] we see that most of the M stars lie at the bottom of 
the izK colour plane. We identify these M stars by assuming their 
typical NIR colour, J - K < 1.3. However, this means that our 
selection criteria must involve J band data and would reduce the 
sky coverage due to the data availability. Therefore we choose to 
exclude these M stars by applying a cut in izK colour plane with 
the condition z- K- 0.9(i - z) > -0.2 in Eq.0 

All magnitudes and colours are given in SDSS AB system 
and are corrected for extinction usi ng the Galactic dust map of 
ISchlegel. Finkbeiner. & DavisI ( Il998h . All colours des cribed below 
refer to the differences in 'model' magnitudes (see iLupton et al.l 
I2OOII for a review on model magnitudes) unless otherwise stated. 

Applying the above selection criteria (Eq. [T) on the SDSS 
DR7, we have two main LRG samples with a total observed area 
(after masking) of a 200deg-. The first sample has 130819 LRGs 
candidates with a sky surface density of a 700deg"- and the second 
one 44543 with a sky density of ~ 240deg"^ The 240 deg"^ LRG 
sample was selected in such a way to check if the redshift distri- 
bution implied by cross-correlations is higher than the 700 deg"^ 
LRG sample. 



the first and second LRG samples, respectively. Using a high num- 
ber density random catalogue helps to ensure the extra shot noise 

is reduced as much as possible. 

We also compute w{6) by using the lHamiltorjh993l) estimator 
which does not depend on any normalisation and is given by. 



DD(e)RR(e) 



(4) 



The Landy-Szalay estimator when used with our samples gives 
negligibly different results to the Hamilton estimator. Note that the 
Landy-Szalay estimator is used throughout this work except in 37.11 
where we used both estimators to test for any possible gradient in 
number density of our samples. 

For the com putation of the cross-correlatio ns in ^and 0we 
use the estimator ( IGuo. Zehavi. & Zhend ( l2oTlh ) : 



WcmsAd) ■ 



DcDsie) - DaRs(e) - DsRD{Q) - RaRs(e) 



(5) 



RaRsid) 

where the subscript G and 5 stands for the contribution in the pairs 
of the quantities that are cross-correlated in each case. 

3.2 Error Estimators 

To determine statistical uncertainties in our methods, we used three 
different methods to estimate the errors on our measurements. 
Firstly, we calculated the error on w{d) by using the Poisson es- 
timate 

1 + w(e) 



a-po. 



y/DIXffj 

Secondly, we used the field-to-field error which is given by 



(6) 



3 THE 2-POINT ANGULAR CORRELATION FUNCTION 
MEASUREMENTS AND ERRORS 

3.1 w(0) Estimators 

The probability of finding a galaxy within a solid angle SCI on the 
celestial plane of the sky at a distan ce 6 from a randomly chosen 
object is given bv(e.g. |Peeble3ll980l) 



dP = n[l +w(e)]Sn, 



(2) 



where n is the mean number of objects per unit solid angle. The 
angular two-point correlation function (2PCF) in our case, actually 
calculates the excess probability of finding a galaxy compared to a 
uniform random point process. 

Different estimators can be used to calculate w(0), so to start 
with w e use the minimum variance estimator from lLandv & Szalavl 
( Il993l) . 



Wls(9) = 1 + 



N,AN,.d - l)\ DD(6) 



N{N- 1) / RR(e) 



RR(e) 



(3) 



where 00(6*) is the number of LRG-LRG pairs, DR(e) and RR(e) 
are the numbers of LRG-random and random-random pairs, respec- 
tively with angular separation 9 summed over the entire survey area. 
N,.^i is the total number of random points, N is the total number of 
LRGs and N,.[i/N is the normalisation factor. For our calculation 
we used two LRG samples (as explained in § 12. U with different 
sky density, thus the density of the random catalogue that we use 
is ~ 20 times and ~ 60 times the number of the real galaxies for 



1 



A'- 



-I 



DRi(e) 
DR(e) 



(7) 



where N is the total number of subfields, Wj(d) is an angular corre- 
lation function estimated from the /th subfield and w(d) is measured 
using the entire field. For this method we divide our main sample 
to 36 subfields of equal size ~ 6deg-. We also reduce the number 
of subfields down to 18 with sizes of ~ \2deg- as we want to test 
how the results could deviate by using different sets of subsamples. 
While Stripe 82 has only ~ 2.5deg height, our subfields with their 
~ 2.5deg and ~ 5deg widths are a reasonable size for estimating 
the correlation function up to scales of 1 - 2deg. 

Our final method is jackknife resampling, which is actu- 
ally a bootstrap method. This technique has been widely used 
in clustering analysis s t udies with correlation functions (see , 
e.g IScrantoii et al.l |2002|: IZehavi et al.l l2005al : IRoss et al.l IIoOtI : 
iNorberg et aLlbOOgl : ISawangwit et al.ll201 ih . The jackknife errors 
are computed using the deviation of the w(6) measured from the 
combined 35 subfields out of the 36 subfields (or 17 out of 18 when 
18 subfields are used). The subfields are the same as used for the es- 
timation of \he field-to-field error above. w(0) is calculated repeat- 
edly, each time leaving out a different subfield and hence results in 
a total 36 (or 18) measurements. The jackknife error is then 



,(e) = £ 



DR{G) 



w,v(£»)- w(0)]-. 



(8) 



where Wji (9) is a measurement using the whole sample except the 
ith subfield and DRi,{e)/DR(e) is approximately 35/36 (or 17/18) 
with slight variation depending on the size of resampling field. A 
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comparison of the error estimators can be seen in Fig. |2] Poisson 
errors are found to be much smaller compared to jackknife errors 
particularly at larger scales. Field-to-field errors give similar results 
as jackknife errors, except aX 6 > 10' where the FtF eiTors under- 
estimate the true eiTor due to missing cross-field pairs. Since the 
jackknife errors are better at a scale of order 100' which are of 
prime interest here, these are the error estimators that will be used 
in this work unless otherwise stated. 

When calculated in small survey areas, w{9) can be affected 
be an 'integral constraint', ic. Normally wiff) has a positive sig- 
nal at small scales and if the surveyed area is suffi ciently small, this 
will c ause a negative bias in w{6) at largest scales jOroth & PeeblesI 
fl977l) . i.e. We^jiO) — wiff) — ic. The inte gral constraint can be calcu- 
lated from (see e.g. lRoche & Ealesll999l) : 



zmff)w,„„M{ff) 



(9) 



where for the w^odeiiff) we assume the standard ACDM model in the 
linear regime 05.31 . No integral constraint is initially applied to our 
full sample results as the expected magnitude of ic is smaller than 
the w{ff) amplitudes at scales analysed in this paper. This position 
will be reviewed when we move on to discuss models with excess 
power at large scales in ^ 

To provide robust and accurate results from the correlation 
functions, we are also interested in model fitting to the observed 
w{9) (see in 35.21 35.4l and 35. 3t . Hence, for model fitting we will 
use the covariance matrix, which is calculated by: 



(10) 



where the w;(0,) is the correlation function measurement value ex- 
cluding the i''' subsample and the factor A' - 1 corrects fr om the 
fact that the realizations are not independent (Myers et al]|2007l; 
Norberg et al.ll2009HRoss. Percival. & Brunneij|2010<:iCrocce et al.l 
201 ll : ISawangwit et al.ll20 1 ij). The jackknife errors are the square- 
root of the diagonal elements of the covariance matrix, so we can 
now calculate the correlation coeflficient, which is defined in terms 
of the covariance. 
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(11) 



where tr^ = C„ (see Fig.[3]l. We can see that the bins are strongly 
comlated at large scales. The covariance matrix is more stable 
when we use 36 Jackknife subfields instead of 18, so we will use 
only the covariance matrix for the case of 36 subfields. 



3.3 Angular Mask and Random Catalogue 

To measure the observed angular correlation function we must 
compare the actual galaxy distribution with a catalogue of ran- 
domly distributed points. The random catalogue must follow the 
same geometry as the real galaxy catalogue, so for this reason 
we apply the same angular mask. The mask is constructed from 
'BEST' DR7 imaging sky coveragfl Furthermore, regions ex- 
cluded in the quality holes defined as 'BLEEDING', 'TRAIL ', 
'BRIGHT.STAR ' and 'HOLE'. The majority of the holes in the 
angular mask is from the lack of K coverage in Stripe 82. The final 
mask is applied to both our data and random catalogue (see Fig.|4ll. 




Poisson/Jackknife 18 subfields 

1 i (1 ,' Jackknife 18 subfields. 
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Figure 2. Comparison of the measured error ratios of the Jackknife, field- 
to-field and the Poisson errors for the w(0) measurements of the 700 deg^- 
Stripe 82 LRG sample. Two diff'erent resampling sets have been used for 
the Jackknife and field-to-field errors, the first one based on 36 subfields 
and the second from 18 subfields. 
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Figure 3. The correlation coefiicients r,j, showing the level of congelation 
between each angular separation bin for the 700 deg"^ Stripe 82 LRG sam- 
ple as calculated by using 36 subfields. 




Figure 4. A fraction of the total ~ 200 deg observed area in Stripe 82. 
LRG candidates (red) and random objects (blue), follow the same angular 
selection. Empty sky patches resulting from the lack of K-band coverage in 
the combined optical-IR data. 
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Figure 5. Noitnalised redshift distributions of MEGAz-LRGs, DEEP2 
galaxies and SDSS QSOs in Stripe 82 tliat are used in tlie cross-correlations 
with the LRG samples. 



For generating the randomly distributed galaxies/points, we 
tried two different ways in order to modulate the surface density 
of the random points to follow the number density and the selec- 
tion function of the real data. The selection function of the random 
catalogue mimics only the angular selection of the real data. 

For the first method, we use a uniform density for the random 
points across the Stripe 82 area, so the normalization factor, N^dlN, 
to be ~ 20 and ~ 60 for the 700 deg - and the 240 deg"^ LRG 
samples, respectively. A second random catalogue was created by 
dividing Stripe 82 into six smaller subfields (15 x 2.5deg^ each) and 
normalizing the density of random points to the density of galaxies 
within each subfield. The difference between the measured angular 
coiTelation function when we use the 'global' or the 'local' random 
catalogue is negligible. We will use the 'gl obal' random catalo gue 
for the clustering analysis. A W-trees code jMoore et al.l200l[) has 
been used to minimise the computation time required in the pair 
counting procedure. 



4 LRG N(Z) VIA CROSS-CORRELATIONS 

Even if the redshift of individual galaxies is not available, the 3-D 
clustering information can yet be recovered if the sample's redshift 
distribution, n(z), is known. This c an be achieved using Limber's 
inversion equation ( lLimbej|l953l) which can project the spatial 
galaxy correlation function, ^(r), to the angular conelation func- 
tion given the n(z) of the sample: 



w(e) = 2 



r dxfix)- r 

Jo Jo 



du ^(r ■■ 



(12) 



where f(x) is the galaxy redshift selection function. For our pho- 
tometric selected LRG samples, only a very small fraction has a 
measured redshift, thus it is vital to estimate the n(z) of the Stripe 
82 LRG samples. 

One method for estimating the redshift distribution of the sam- 
ple could be based on the various popular programs that derive pho- 
tometric redshifts (photo-z's). Photo-z estimates are based on the 
deep multi-band photometry coverage, and work by tracing some 
specific spectral features across the combination of filters which 
are then compared with difl"erent type of objects SED templates. 



Indeed, our izK selection is a rough photo-z cut as we follow the 
movement of the 4000A break across the selected bands. In or- 
der to use the angular correlation function and the information that 
is encoded we need the n(z) of our sample, hence we follow the 
technique of iNewmanI j2008l) for reconstructing the LRG redshift 
distribution from cross-correlations. 



4.1 Redshift distribution reconstruction 

We employ Newman's method, which is about determining 
the underlying redshift distribution of a sample of objects (LRGs 
in our case) through cross-correlation with a sample of known red- 
shift distribution. By cross-correlating the sample (or samples) with 
known redshift and the sample under consideration, if both samples 
lie at the same distance, this will give a strong clustering signal. If 
the two samples that we are cross-correlating are separated and are 
at different z distances, no cross-coixelation signal will result. Thus, 
through the cross-correlations we can infer our photometrically se- 
lected LRG sa mple z ranges. 

Following lNewmaiT ( l2008l) the probability distribution func- 
tion of the redshift of the Stripe 82 LRG samples, (fip(z), is: 



^p(z) = w{z) 



3-7 dA(zf-dI/dz 



2n 



(13) 



where w{z) is the integrated cross correlation function, w,p(0, z), 
of the LRG photometric samples with the samples of known spec- 
troscopic redshift (see ^A?2S . H(y) = r(l/2)r((7- l)/2)/r(7/2) 
where T(x) is the Gamma function, is the comoving angular dis- 
tance and dl is the comoving distance at redshift z. The comoving 
distance r,„„j. corresponds to the maximum angle at given redshift, 
which must be large enough to avoid nonlinear biasing effects. 

To derive ippiz) via Eq. [T3] we must estimate Wsp(9,z) ~ 
(ppiz) rl^'^sp, since the angular size distance, d^iz) and the comoving 
distance l{z) are given by the assumed cosmology. Thus we now 
require only knowledge of the 7,,, and ro ,,, parameters as function 
of redshift. Fortunately under the assumption of linear biasing, the 
cross-correlation of the two samples under consideration is the re- 
sult of the geometric mean of the autocorrelation functions of the 
samples, i.e. = (f.,,5<f,,p)i , hence we can use the information pro- 
vided by autocorrelation measurements for each sample to break 
the degeneracy between correlation strength and redshift distribu- 
tion. 

Newman investigates the effect of systematics such as: dif- 
ferent cosmologies, bias evolution, errors from the autocorrelation 
measurements and field-to-field zero points variations in the final 
redshift probability distribution result. These issues could be more 
important in the case of future photometric surveys aimed at plac- 
ing constraints on the equation of dark energy. 



4.2 Cross-Correlation data sets 

Newman's angular cross-correlation technique requires the use of a 
data sample with known spectroscopic, or sufficiently accurate pho- 
tometric, redshifts. For this reason we use a variety of samples with 
confirmed spectroscopic and photometric redshifts for the cross- 
coiTelations with Stripe 8 2 LRGs. The data samp les that we use 
are: DEEP2 DR3 gal axies l lDavis et al.l200ll2007l) . MegaZ-LRGs 



ega 

jCollister et al.l2007h. SDSS DR6 QSOs ^Richards et al.l2009l) and 

^ 120101) . In 



SDSS DR7 QSOs l lSchneider et al 



Fig. [5] we show the 
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Table 1. Number of objects in each separate redshift-bin used for the cross-correlations with Stripe 82 LRGs 



redshift 


DEEP2 


sample 

MegaZ-LRGs DR6 Photometric Sample 


DR7 Spectroscopic sample 


0.4 - 0.6 




30503 436 


456 


0.6 - 0.8 


3152 


695 


526 


0.8- 1.0 


5512 


1199 


547 


1.0- 1.2 


3620 


1630 


729 


1.2- 1.4 




1312 


820 


1.4- 1.6 




2646 


854 


1.6-1.8 




1193 


803 


1.8-2.0 




1990 


668 





Figure 7. Cross-coiTelation measurements of the 240 deg ^ and 700 deg ~ Stripe 82 LRG samples with DEEP2 galaxies in (a) and (b), respectively. Uncer- 
tainties are Icrjackknife errors. 
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Figure 8. Cross-correlation measurements of the 700 deg"~ (green dia- 
mond) and 240 deg"' (purple star) Stripe 82 LRGs with MEGAz-LRGs. 
along with lo" jackknife errors. 



4.2.2 DEEP2 Sample 

The next sample of galax ies that we use is DEEP2 DR3 galaxies 
l lDavis etalj[200l |2007|) . The survey coverage in Stripe 82 is 
~ 1.7 deg^ with i < 24. Galaxies in DEEP2 are split in three 
redshift bins with 0.2 step in the redshift range 0.6 ^ z < 1.2. The 
redshift distribution of the DEEP2 DR3 sample is shown in Fig. [5] 
with 12,284 galaxies in total. In Fig. [7] we show the results of the 
cross-correlations of the 700 deg"^ and 240 deg"^ LRG samples 
with the DEEP2 galaxies in the three aforementioned redshift bins. 



4.2.3 MegaZ-LRG sample 

The last sample that we use are LRGs fro m the MegaZ-LRG pho- 
tometric catalogue ( ICollisteretai1l2007l) . MegaZ-LRGs are used 
only in the redshift range of 0.4 ^ z ^ 0.6 with / < 20. This sam- 
ple offers us the ability to check the clustering properties of our 
high-redshift LRG candidates with another sample of LRGs. The 
total number of MEGAz LRGs that we use for cross-correlations 
is 30,503. In Fig. [8] are shown the cross-correlations between the 
Stripe 82 LRGs and the MEGAz LRGs. 



normalised redshift distributions of all the samples and in Table [T] 
we present the number of objects in each redshift bin. 

By using the above data sets for cross-correlation we satisfy 
the principal requirements of Newman's method, with the most im- 
portant being that the sky coverage of the data sets overlap the 
Stripe 82 LRGs. It must be mentioned though that not all the red- 
shift surveys have the same sky coverage as Stripe 82 LRGs, so we 
reconstruct two redshift distributions via the cross-correlations pro- 
viding us with the opportunity to check how much the n(z) cross- 
correlation technique is affected by area selection. One n{z) is re- 
constructed by using all the data sets, the other n{z) by using only 
SDSS QSOs in the cross-correlations. 



4.2.1 SDSS DR6 & DR7 QSOs 

QSO surveys are the main samples that we used for our cross- 
correlation measurements and they span the redshift range 0.4 < 
z < 2.0. When we refer to QSO data sets, we separate them into 
spectroscopic and photometric samples. 

For the spectroscopic QSO sample we use the fifth edition 
of the SDSS Quasar C atalog, which is based on the SDSS DR7 
jSchneider etaU l2010h . The original data set contains 105,783 
spectroscopically confirmed QSOs, from which only 5,403 in 
Stripe 82 have been used at 0.4 < z < 2.0 for cross-coiTelations 
(Table[3 with / < 22 (~ 28% of QSOs at i > 20). 

The photometric QSO sample comes from the photometric 
imaging data of the SDSS DR6 (Richards et al, 200.^. The parent 
catalogue contains ~ 1, 000, 000 QSOs candidates from which we 
use 11,101 with i < 21.3 in Stripe 82 and in the same redshift range 
as the spectroscopic QSOs. 

In Fig. [6] we plot the cross-correlations between the Stripe 
82 LRGs and the SDSS QSOs. We show only the case for cross- 
correlations of the 700 deg"' Stripe 82 LRG sample with the spec- 
troscopic and photometric SDSS QSOs. Cross-correlation with the 
240 deg"* LRG sample does not differ much. Errors shown here 
and for the other cross-correlation cases are jackknife errors. 



4.3 Cross-Correlation results for n(z) 

Having estimated the clustering signal from the cross-correlations 
of the above samples, we proceed to the reconstruction of the red- 
shift distribution of the photometrically selected Stripe 82 LRG 
candidates. To estimate the probability distribution function of the 
redshift, (f>p(z), for the high-z LRG candidates we use equation il3l . 
The pair-weighted clustering signal of the cross-correlations has 
been integrated up to a 6' for each redshift bin. 

In Fig. |9] we can see the two cases of the estimated proba- 
bility distribution function of the redshift for the high-z LRG can- 
didates. For the first case, (f>p(z) has been estimated by using the 
spectroscopic SDSS QSOs whereas in the other case, (^,,(z) is es- 
timated using only the photometric SDSS QSOs (DEEP2 galaxies 
and MEGAz-LRGs are also always used). For both cases we plot 
the en'ors estimated for each point in the redshift bin from the con- 
tributed cross-coiTelated sample. 

To estimate the redshift distribution, n(z), we use the weighted 
mean for the (f),,(z) in each redshift bin, calculated through : 

k 

n(z) = , (14) 

2(1/-?) 

1=1 

where k is the total number of bins at that redshift, (pp(i) is the mea- 
sured probability distribution function of each cross-correlation 
data set in the ith bin and cr, the error on that measurement. 

The spectroscopic QSO (/>p(z) in Fig. compared to the 
photo-z case in Fig.|9]5, gives increased probability at z ~ 1. This 
may be explained by the SDSS QSO spectroscopic redshifts being 
more precise. For this reason, in our analysis and in fitting mod- 
els to our w{d) results, we will use only the spectroscopic n(z) for 
higher accuracy. 

In Fig. [To] we plot the normalized redshift distribution of the 
240 deg"^ and 700 deg"^ LRGs samples as calculated from Eq. [T3] 
- 1141 When we selected the two LRG samples from the izK colour- 
plane, we applied a redder selection for the 240 deg"^ sample (see 
Eq. \Q, aiming for a sample with a slightly higher redshift peak in 
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Figure 9. (a) The probability distribution function of the redshift, ifip(z), of the 700 deg"- and 240 deg"- Stripe 82 LRGs as estimated through cross-correlations 
with MEGAz-LRGs, DEEP2 galaxies and spectroscopic SDSS QSOs. (b) Same as in (a) but now using photometric SDSS QSOs instead of spectroscopic in 
the cross-correlations. Error bars shown in both cases are lo" jackknife summed up to 6'. 
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Figure 10. Weighted normalised redshift distribution of the Stripe 82 LRGs 
candidate samples when we use the spectroscopic SDSS QSOs along with 
the DEEP2 and MEGAz-LRG data sets. As expected the 700 deg"^ sample 
(solid green line) n(z) peak is lower when compared with the 240 deg"^ 
sample (dashed blue Une). 



the distribution as predicted from the evolutionary tracks in Fig.[T] 
This small difference may be seen between the spectroscopic n(z) 
of the 700 deg"^ and 240 deg"- samples where the bluer cut has 
an average of z ~ 1 where for the redder sample the average is 
z~ 1.1. But since the 700 deg"^ LRG sample has higher statistical 
accuracy in the n(z) determination, the majority of our analysis will 
be focused in this sample. 



5 RESULTS 

5.1 Measured w{ff) and comparisons 

In Fig. [TT] we compare the observed ang ular correlation function 
of the 700 deg-2 LRG in Stripe 82 with ISawangwit et all l [20Tlh 
results. The w{9) measurements are presented with Icr Jackknife 
errors. 



10.000 E 



1.000 E 



Stripe82 LRGs 
AAn LRGs □ 
2SL.'\0 Ll^G 
SDSS LRGs o 



0.001 L 
0.1 




1.0 10.0 
G/ acmin 



100.0 



Figure 11. The angular correlation function, w{6), from the 700 deg ^ 
Stripe 82 LRGs (star), AAQ. LRGs (square), 2SLAQ LRGs (triangle) and 
SDSS LRGs (diamond). At small scales all of the measurements show sim- 
ilar clustering behaviour, but at large scales the Stripe 82 clustering slope 
appears to be flatter than the lower z. samples. 



The work of lSawangwit et all involved three LRG data sets at 
z< 1 : 

(i) SDSS LRGs at z ~ 0.35 

(ii) 2SLAQ LRGs z ~ 0.55 

(iii) AAn LRGs z ~ 0.68 

From Fig.[TT]we can see that at small scales, 9 < 1', the clus- 
tering trend for all the samples is similar but with decreasing ampli- 
tude for increasing redshift. At larger scales, we note that the w(9) 
of the Stripe 82 LRGs seems to have a flatter slope than the other 
samples, departing from the expected behaviour for the correlation 
function. 

Further com parisons below with the LRG clustering results of 
ISawangwit et al.l will focus on the slope and amplitude of the w(6) 



Clustering analysis of high-redshift Luminous Red Galaxies in Stripe 82 9 



results, with an initial view to interpret any changes in terms of 
evolution. It is therefore of interest to see how the Stripe 82 sample 
match to the LRG samples used in previous studies in terms of 
luminosity and comoving space density. 

A pair-weighted galaxy number density is given by (see e.g. 
IRoss & Brunnerll2009l) : 



H(zMz) 



-Q„t,,cP(z) "" ' ^^^^ 
where n„iys is the observed area of the sky, l(z) is the comoving 
distance to redshift z and c is the speed of light. The observed space 
density for the VOOdeg"- Stripe 82 sample is found to be » 3.20 ± 
0.16 X lO^-^h'Mpc-^ The quoted la error has been estimated from 
the difference of the number density as calculated through Eq. [TS] 
and by converting Fig.llOlinto a plot of number density as a function 
of z (by dividing its bin by its corresponding volume). 

Within the uncertainties of our n(z), the 700deg"^ sample ap- 
pears to have similar space density to that of the AASl LRG sam- 
ple (see Table |2] in 35. 2t . However, in this study we do not yet 
have redshift information for individual LRGs, not even for a sub- 
set of the sample. Hence it is more uncertain i f our sample has 
simila r luminosity as the LRG samples used bv ISawangwit et al] 
( I2OIII) . We therefore take the fact that the samples are number- 
density matched to imply that they are also approximately luminos- 
ity matched which may tur n out to be a reasonable assumption (see 
e.g. lSawangwit et alj201ll) . This then should enable us to compare 
the clustering slopes and amplitudes of the AACl and Stripe 82 and 
infer any evolution independently of luminosity dependence. 

5.2 w{9) and power-law fits 

Our first aim here is to fit power-laws to the Stripe 82 w(9) to pro- 
vide a simple parameterisation of the results. Our second aim is to 
make comparisons of the 3-D correlation amplitudes and slopes to 
measure evolution. Both aims will require application of Limber's 
formula to relate the 2-D and 3-D correlation functions. 

We begin by noting that the simplest function fitted to corre- 
lation functions is a single power-law with amplitude rg and slope 
7. In previous studies, the spatial coiTelation function has been fre- 
quently described by a power-law of the form: 



(16) 



The angular correlation function as a projection of ^(r) can be writ- 
ten as w(d) = Q'S'"'', commonly with a slope fixed at 7 = 1.8. The 
amplitude of the angular correlation function, a, can be related with 
the correlatio n length rp through Limber's fo rmula (Eq. I12t using 
the equation ( jBlake. Collister. & Laiiavll2008l) : 



C P 



dz n(z) 



xizY 



(17) 



where n(z) is the redshift distribution, x(z) is the comoving ra- 
dial coordinate at redshift z and the numerical factor Cy = 

ra)r(f-i)/r(^). 

A deviation from a single power l aw at ~ Ih 'Mpc 
has been measured in previous studies (Shanks etal. ' 19831; 



Blake. Collister. & Lahavl .2008: Ross et al. 2008; Kim et al. 2011; 
Sawangwit et al.l 201 ll) and can be explained by the the 1-halo 
and 2-halo terms imprinted in the clustering signal under the as- 
sumption of the halo model (see 35.4b . To parameterise the clus- 
tering characteristics of our sample, we fit a single-power law and 



1.00 



£ 0.10 



0.01 



single power-law 
double power-law 

w(e) + loJK error 



• — I- 



-i — i — I 



Figure 12. The best-fit single power law (diamond) and double power law 
(triangle), for the 700 deg"- LRGs candidates overplotted on the angu- 
lar congelation function (square) with the 1 <t Jackknife error. Lower panel 
shows the fitting residuals. 



a double-power law to our measured angular correlation function. 
The double power-law form is given as: 



W2(e) 



eo.i 
9 



1-71 



(6 < Ot) 



l-yn 



(9 > Ob) 



(18) 



(19) 



with Ob to be the break point at a 1.2' where the power-law slope 
changes from being steeper at small scales (< 1.2'), to flatter at 
large scales. 

The power-laws are fitted in the range OA' < 6 < 30' using the 
-minimization with the full covariance matrix constructed from 
the jackknife resampling (see 33. 2t : 



(20) 



i.j=i 



where A' is the number of angular bins, Aw(dj) is the difference 
between the measured angular correlation function and the model 
for the ith bin, and Cjj is the inverse of the covariance matrix. 

For the single power-law, our best-fit spatial clustering length 
and clustering slope pair from Limber's formula are measured to 
be ro = 7.54 ± 0.16h"'Mpc and 7 = 2.01 ± 0.01 with associated 
reduced = 5.89. The ro - 7 pairs for the double power-law are 
ro.i = 7.63 ± 0.27h"'Mpc and y, = 2.01 ± 0.02 at small scales and 
ro,2 = 9.92 ± 0.40h"'Mpc and 72 = 1.64 ± 0.04 at large scales with 
a reduced = 3.65. From the intersection of the 2 power law for 
f(r), we have calculated the break scale, rb = 2.38h"'Mpc. This is 
higher than the rb = 1.3 - 2.2h ~'Mpc estimated from the SDSS, 
2SLAQ and AAQ LRG surveys jSawangwit et alj201 ll) . 

In Fig. [T2] we show the data points including the lo" Jack- 
knife en'ors with the best-fitting power laws where the largest 
scale considered in the fitting was < 30', which corresponds to 
r < 20h-^Mpc at z ~ 1 for the 700 deg-^ LRG sample. Fig. [HI 
confirms that the double power-law clearly gives a better fit to the 
data than the single power-law. Note that in the case of the single 
power-law and the double power-law at small scales, our results 



10 A^. Nikoloudakis et al. 

Table 2. Best fit parameters for the single and double power-law fits to the angular correlation function. 



Sample 


z 


"g 


Sin{ 


jle power-law 






Double power-law 












r 


ro(/i"' Mpc) 




ri,2 


'•o.i.2(/i"' Mpc) 


/■(,(/) ' Mpc) 


^ red 


AAn LRGs 


0.68 


2.7 X 10-* 


1.96 ±0.01 


7.56 ± 0.03 


42.8 


2.14 ±0.01 


5.96 ± 0.03 


1.3 


3.4 














1.81 ±0.02 


7.84 ± 0.04 






Stripe82 LRGs 


1.0 


3.20 ± 0.16 X \0-^ 


2.01 ± 0.01 


7.54 ±0.16 


5.89 


2.01 ± 0.02 


7.63 ± 0.27 


2.38 


3.65 


{100 deg-^ ) 












1.64 ±0.04 


9.92 ± 0.40 








0.1 1.0 10.0 100.0 

O/arcmin 



Figure 13. a) The AAn LRG raw w{8) measurements (red square) with pre- 
dictions from comoving evolution model (dashed red line), using the best-fit 
double power-law tq- y values with Limber's formula as lSawangwit et aU 
ilOl ll) calculated. We then evolve the AAfi best-fits utilising the estimated 
700 deg"^ Stripe 82 LRG n(z) under the assumption of comoving evolu- 
tion (dashed-dot black line) clustering. The observed Stripe 82 LRG w(8) 
is shown as well (blue star), b) Same raw measurements as above, but 
now compared to the virialised evolution clustering model. Stripe 82 LRG 
w{6) measurements are described more accurate with comoving evolution 
at small and large scales compared to virialised evolution as it can been seen 
from the lower panel, where are plotted the residuals of the observed Stripe 
82 w(ff) versus the comoving evolution (black stai') and virialised evolution 
(green star) models, respectively. 



give ?■() - 7 values consistent with outcomes from previous studies. 
However, at large scales the Stripe 82 slope (72 = 1-64 ± 0.04) is 
significantly flatter than the AAfi result (72 = 1.81 ± 0.02). 

Fig. [T3] shows the double power -law fits for AAQ (dashed 
red lines) taken from lSawangwit et al.l and then evolved (black and 
green dot-dashed lines) to the Stripe 82 depth using Eq. [17] un- 
der the assumptions of comoving and virialised clustering, respec- 
tively. We shall interpret the amplitude scaling in the discussion of 
evolution in 36.1 H ater. At this point we again note that the biggest 
discrepancy seems to be at large scales where the Stripe 82 slope is 
increasingly too flat relative to the AAfl result. Fitted parameters 
are given in Table |2l where the best-fit power-law parameters for 



the AAfi LRG sample jSawangwit et alj[201 ll) are also presented 

for comparison. 

We note here that iKim et al] j201 ih studied the clustering 
of extreme red objects (EROs) at 1 < z < 2 in the SA22 
field and they report a sim ilar change of the large scale slope. 
iGonzalez-Perez et al. I( |20T l|) tri ed to fit clu stering predictions from 
semi-analytic simulations to the lKim et alj FRO w{6) but found that 
the model undeipredicts the clustering at large scales. 

5.3 ACDM model fitting in tlie linear regime 

Since the standard ACDM model w as found to give a good fit to 
the lower redshift LRG samples of ISawangwit et al.l ( l201lb . we 
now check to see whether the flatter large-scale slope of the Stripe 
82 LRG w(6) leads to a statistically significant discrepancy with 
the ACDM model aXz ~ 1 . We generate matter power spectra us- 
ing the 'CAMB' software l lLewis. Challinor. & LasenbvlF2000l) . in- 
cluding the case of non-linear growth of str ucture correction. For 
this reason we use the 'HALOFIT' routine jSmith et aLll2003h in 
'CAMB'. Our models assume a ACDM Universe with Ha = 0.73, 
Q.,„ = 0.27, /4„,.™„ = 0.167, CT-g = 0.8, h = 0.7 and n., = 0.95. 
Then we transform the matter power spectra to obtain the matter 
correlation function, <fm(r), using: 

ryiCO Slflkf 

^nM=— p,„(k)e—^dk. (21) 

2?r2 Jo kr 

The relationship between the galaxy clustering and the under- 
lying dark-matter clustering is given by the bias, bg : 



As we are interested in the linear regime, we fit the projected 
^n,(r) to the Stripe 82 LRG w(e) in the range 4' < < 45', corre- 
sponding to comoving separations 3 < r <, 30h"'Mpc. By fitting 
the model predictions to the measured w{6) it will result with the 
best linear bias factor, the only free parameter in this case . For our 
fitting, the ^^-minimization with the full covariance matrix con- 
structed from the jackknife resampling (see 33. 2t has been used. 

The best-fit linear bias parameter is estimated to be = 
2.74±0.07 withx^.^j = 5.09. The upper limit of our fitted range in 9 
was varied, while the lower limit stayed constant to avoid any con- 
tribution from the non-linear regime. Thus, for the range ~ 4' - 30' 
the best-fit bias is b = 2.8 ± 0.08 with = 4.72 and at ~ 4' - 60' 
is b = 2.69 ± 0.07 with xia = 5-18. In Fig.[T4]we plot the LRG 
w(d) with the Icr error and the ACDM model with the best-fit bias. 
For low values of the upper limit of the fitting range, the measured 
biases are in approximate agreement with other results in the lit- 
erature. But in terms of the flat slope of w(6) at large scales, the 
standard ACDM linear model is inconsistent with the data at the 
2 - 3iT level. One of the aims of the next section will be to see if 
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Figure 14. The best-fit spatially flat, ACDM model assuming n,„ = 0.27 
compared to the observed n'(6) of Stripe 82 700 deg"- LRGs in the linear 
regime. The standard model cannot explain the large scale power excess in 
the angular correlation function of the Stripe 82 LRGs. The shaded area 
corresponds to ±1it jackknife error Also shown is a spatially flat ACDM 
model with the same parameters as before except for a lower value of £1,,, = 
0.2 and an arbitrary normalisation. The VL,„ = 0.2 model appears to give a 
better fit than the standard fl,,, = 0.27 model. 



a HOD model can explain the flat large-scale slope of the z ~ I 
Stripe 82 LRGs. 



5.4 Halo model analysis 

We are going to us e the approach of the halo model (see 
ICoorav & Sh"eflill2002l for a review) of galaxy clustering to finally 
fit our angular correlation function results. Under the halo-model 
framework we can examine the way the dark matter haloes are 
populated by galaxies through the Halo Occupation Distribution 
(HOD). Various studies have used this model to fit their results (e.g. 
Masiedi et al.' '200f); 'White et al."2007'; 'Blake, Collister, & Lahavl 
2008; Wake et al. 2008; Brown et al. 2008; Ross & Brunner 2009; 



Zheng et al.l |2009| ; ISawangwit et al.l l201ll ; I Gonzalez-Perez et al.l 
201 ih as a way to explain the galaxy correlation function and 



gain insight into their evolution. Specifically, we shall investigate 
whether the HOD model may be able to explain the flatter slope of 
the correlation function observed here. 

In the halo model, the clustering of galaxies is expressed by 
the contribution of number of pairs of galaxies within the same 
dark matter halo (one-halo term, ^1) and to pairs of galaxies in two 
separate haloes (two-halo term) : 



(23) 



The 1-halo term dominates on small scales < IMpc. 

The fundamental ingridient in the HOD formalism of galaxy 
bias is the probability distribution P{N\M), for the number of galax- 
ies A' to hosted by a dark matter halo as a function of its mass M. 

We use the so-called centre-satel l ite three-parameter HOD 
model (e .g. 'Seo. Eisens tein. & Zehavil |2008| ; IWake et all |2008| ; 
ISawangwi t et al. 201 1) which distinguishes between the central 
galaxy and the sa tellites in a halo. This separation has been shown 
in simulatations jKravtsov et"!!] |2004|) and has been commonly 
used in semi -analytic galaxy formation models in the last years 
( lBaughll2006l) . 



Different HODs are applied for the central and satellite galax- 
ies. We assume that only haloes which host a central galaxy are 
able to host satellite galaxies. The fraction of haloes of mass M 
with centrals is modelled as: 



{N^\M) = exp 



/' -M„i„ \ 



(24) 



In such halo es, the number of sa tellite galaxies follows a Poisson 
distribution dKravtsov et alj2004l) with mean: 



{NAM)) 



(25) 



To describe the distribution of the satellite galaxies around the hal o 
centre we use the NFW profile ( INavarro. Frenk. & Whitelll997b . 
So, the mean number of galaxies residing in a halo of mass M is: 

(N\M) = (N,\M) X (1 + (N,\M) . (26) 

and the predicted galaxy number density from the HOD is then: 

Kg = J dMn(M){N\M) (27) 

where n( M) is the halo mass func tion, where in our case we use the 
model of lSheth & LemsorJ l ll999l) . 

From the HOD we can derive useful quantities which are the 
central fraction : 

r dMn(M) <Nc(M)) 

f = d 

dM n(M) <Nc(M)) [1 + (K(M))] ' 

and the satellite fraction of the galaxy population: 



(28) 



(M){N,(M)) W|M), 



(29) 



as Fot, = 1 - F^.g„. We can also determine the effective mass, M^ff, 
of the HOD: 



effect 
■ J dM m 



Meff = — I dM n{M)M {N\M) , 

and the effective large-scale bias: 
I 



{M)h{M) {N\M) , 



(30) 



(31) 



where b{M) is the halo bias, for which we u se the ellipsoidal col- 
lapse mode l of Sheth, Mo, & TormenI 1 I2OOII) and the improved pa- 
rameters of iTinker et al. ( 20051) . 

As the galaxy congelation function is the Fourier transform of 
the power spectrum, the 1-halo term and the 2-halo term of the 
clustering functions can be written as: 



P{k) = P,^{k) + P.Ak). 



(32) 



Moreover the 1-halo term can be distinguished from the contribu- 
tion of the cent ral-satellite pairs, P^sjk ), and satellite-satellite pairs, 
P,M), (see e.g. lSkibba & Shethll2009l) : 



Pc.(k) 
and 
P.Ak) ■■ 



fdMm 



(M)2 {N^M) {N,\M) u {k\M) , 



dMn(M) {N^M) (N,\M)- u{k\M)' 



(33) 



(34) 



where M(k|M) is the NFW density profile in Fourier space 
and we have simplified the number of satellite-satellite pairs 
{NsiNs - 1)|M) to {Ns\M)^ since the satellites are Poisson- 
distributed. 
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Table 3. Best-fit HOD parameters. 



Sample z Mi a rig M^ff F,.^t hm X\ 

(lO'^/i-'Mo) (lO'3/i-lMo) (W-^lr'Mpc-^) (IO'^/j-'Mq) (percent) 





0.68 


1.02 


± 


0.03 


12.6 


± 


1.0 


1.50 


± 


0.03 


Stripe82 (10') 


1.0 


3.09 


± 


0.75 


30.2 


± 


6.7 


2.38 


± 


0.12 


Stripe82 (30') 


1.0 


2.57 


± 


0.31 


25.7 


± 


3.1 


2.28 


± 


0.04 


Stripe82 (45') 


1.0 


2.19 


± 


0.63 


21.9 


± 


5.6 


2.24 


± 


0.12 


Stripe82 (60') 


1.0 


2.19 


± 


0.21 


21.9 


± 


2.1 


2.25 


± 


0.05 



3.1 ±0.4 


3.0 ±0.1 


9.0 ±0.09 


2.08 ± 0.03 


13.6 


0.5 ± 0.3 


4.0 ± 0.6 


2.13 ± 1.0 


3.01 ±0.21 


2.4 


0.6 ± 0.2 


3.6 ±0.5 


2.62 ± 0.07 


2.90 ±0.15 


2.3 


0.8 ± 0.3 


3.3 ±0.6 


3.17 ±0.10 


2.81 ±0.18 


3.1 


0.8 ± 0.2 


3.3 ± 0.3 


3.17 ±0.08 


2.81 ±0.10 


3.6 




(a) (b) 

Figure 15. (a) The mean number of LRGs per halo as a function of halo mass at j = 1 . The total, central and sateUite contiibutions are shown by the solid, 
dashed and dotted lines, respectively, (b) The measured angular coiTelation function w(6) for the 700 deg"- LRG sample with the best HOD fit (black star). 
The 1, 2 and 3(T Jackknife errors are shown in red, blue and green, respectively. 



The 2-halo term is evaluated as: 
1 



P^(k) X 



Jo 



Mum('-) 



dM n(M)b(M, r) {N(M)) u(k, M) 



(35) 



where Pmih) is a non-linear matter power spectrum. We derive 
the mass lirtiit, Mii ,n(r), using the '«^-matched' approximation of 
jTinker et al.]|2005l) , which accounts the effect of halo exclusion; 
different haloes can not overlap, n'^ is th e restricted galaxy number 
density (Eq. B13 of lTinker et all j^OOSh ). 

For the scale-dependen t halo bias, b(M, r), we use the model 
given bv lTinker et"alj ( l2005h : 



Z>^(M,.) = .^(M)Ii±HZ^ 
[l+0.69<r„,(r)]' 



(36) 



where fm('') is the non-linear matter correlation function. For the 
2-halo term, we need to correct the galaxy pairs from the restricted 
galaxy density to the entire galaxy population. 

By using Limber's formula to project the predicted spatial 
galaxy correlation function ^(r) to the angular correlation func- 
tion w(9) and we fit for a variety of the three-parameter halo model 
(M„„„, Ml, ff). 

The best-fit model for each of our sample is then determined 
from the minimum value of the^^-statistic using the full covariance 



matrix. We use the full covariance matrix over the range 0.25' < 
6 < 60' in our fitting. Smaller scales are excluded in the fitting 
because any uncertainty in the ^(r) model can have a strong effect 
on w(6) due to the projection. To determine the Icr error on the 
fits, the region of parameter space from the best fits with dx^ < 1 
(Icr for 1 degree of freedom) is considered. For b\i„, M^if, Fsa, and 
rig which depend on all the three main parameters, the considered 
region of the parameter space becomes 6x' < 3.53. 

Fig.ll5b shows the resulting best-fit HOD of the mean number 
of LRGs per halo along with the central and satellite contributions. 
The best-fitting values for M^i„, Mi and a where M„,„ = 2.19 ± 
0.63X lO'^h-'Mo, Ml = 21.9±5.6x 10"h-'Mo and a = 2.24+0.12, 
respectively. The associated values for bn^, M^ff, F^at and «g are 
given in Table [3] 

We see that the {N\M) of the LRGs flatten at unity, as ex- 
pected from the assumption satellite galaxies are hosted by halos 
with central galaxies. The LRGs as expected populate massive dark 
matter haloes with the masses « 10'^^ - 10'''h"'Mo. With the frac- 
tion of LRGs that are satellites being less than 5%, we therefore 
find that > 95% of LRGs are central galaxies in their dark matter 
halo es. The best fit linear bias , b^^ x 2.8, agrees with the prediction 
from lSawangwit et al.l ( 1201 ih in the case of a long lived model for 
the LRGs and indicates that the LRGs are highly biased tracers of 
the clustering pattern. The effective mass, M^g a 3 x lO'^h^'Mg, 
confirms that LRGs are hosted by the most massive dark matter 
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haloes. Despite the fact that we use a higher redshift LRG sample, 
our best-fit HOD parameters are statistically not too dissimilar to 
those found in previous LRG studies (eg see Table [3j- 

In Fig.lISb we show the best-fit model for w(d), compared to 
the data. The first thing we notice is that while at small scales the 
best-fit HOD are in good agreement with the w(d) measurements, 
at large scales the model fits only at 2- 3(T. The flatter slope at large 
scales is responsible for that and we still are not able to say if this 
can be explained by evolution in the linear regime or any kind of 
systematic eifect. In ^we will check systematic errors that could 
affect our results. 

Moreover, due to the high value of the best-fit reduced x~ = 
3.1, we also try to fit the HOD models at different scales by using 
4 different maximum 9 bins of the covariance matrix in our fits, 
which we present in Table |3] The fits at large scales did not im- 
prove and above 45' there was not any change in the best-fit HOD 
measurements. 

Considering the two-halo term in the HOD model, one can 
see that the bias in this regime is mostly scale-independent and the 
correction factor is in fact having the opposite effect on the slope. 
The scale-independent bias is simply the average of the halo bias, 
b(M), weighted by the halo mass function and the mean number 
of galaxies hosted by the corresponding halo. One way to boost 
the large-scale amplitude is to increase M^un and therefore increase 
the mass range of the halo where most galaxies occupy and hence 
linear bias and amplitude of the two-halo term. However, to com- 
pensate for the increase numbers of satellite galaxies (and conse- 
quently small-scale clustering amplitude) one must also increase 
Ml, the mass at which a halo hosts one satellite galaxy on average. 
And in order to produce the overall flatter slope one needs to in- 
crease Mi/M„,i„. However, this would still overpredict the cluster- 
ing amplitude in the intermediate scales, r ~ 5 - 10 h^'Mpc. Note 
that our best-fit HOD gives Mi/M,^i„ a 10, consistent with pre- 
vious results for lower redshift LRGs of dSawangwit et al.lbol ih 
and ( Wake et alj|2008l) . However, as noted earlier including w(9) 
bins at larger and larger scales does not change the best-fit parame- 
ters which means that M[ /M,nin also remains unchanged due to the 
reason discussed above. We therefore conclude that the HOD pre- 
scription in the framework of standard ACDM cannot explain the 
observed large-scale slope in w(6) of the z ~ I LRG sample. 
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linear theory 
stable clustering 
long-lived model 
no-evolution 



<M,> - -22.7 
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Figure 16. The LRG ^20 measurements as a function of redshift and lu- 
minosity from Sawangwit et al. (2011). Lowest redshift data are eai'ly-type 
galaxies from 2dFGRS (Norberg et al. 2002). Stars represent the brighter 
samples (SDSS, 2SLAQ* and AAn*-LRG), where the lower luminosity 
samples, triangles, have been lowered by 0.2 for clarity. The 700 deg"^ 
Stripe 82 LRGs ^20 measurement is at i = 1 (square). 



However, the flat slope beyond Ih^'Mpc causes a highly sig- 
nificant, « 50%, rise in ^20 above the AAQ ^20 as we can see in 
Fig- EH (see also Figs.|13b.b). If we assume that the 2 samples are 
m atched then we would co nclude that all of the models discussed 
bv lSawangwit et aP (I2OI ih were rejected. 

One possibility is that the 700deg"- LRG sample is closer to 
the SDSS and AAQ* LRG space density of 1.1 x lO""* h-^Mpc"' 
because the LRG ^20 fits the extrapolated models better there. If 
so, then this would imply that the Stripe 82 LRG n(z) width was 
underestimated in the cross-correlation procedure and this would 
then increase the deprojected amplitude of ^(r), suggesting that this 
explanation may not work. Similarly a larger correction for stellar 
contamination would also produce a higher Stripe 82 clustering am- 
plitude. We do not believe that loo king further into the evolution of 
the bias jPapageorgiou et al.l2012h and DMH is warranted until we 
understand the flat slope of the Stripe 82 w(d) at large scales. 



6 CLUSTERING EVOLUTION 
6.1 Intermediate scales 

First, we compare the clustering of the z ~ I Stripe 82 LRG 
sample to the lower redshift z ~ 0.68 AACl LRG sample. We re- 
call that these LRG samples have approximately the same space 
density and so should be approximately comparable. We follow 
ISawangwit etal.l ( l201lh andbv using our best-fit ro and 7 we make 
comparison with their data and models via the integrated correla- 
tion function in a 20 h^'Mpc sphere, ^20- 

AA A LRG res ults are described better with the long-lived 
model of lFrvl l ll996h . Fry's model assumes no merging in the clus- 
tering evolution of the galaxies while they move within the gravi- 
tational potential, hence the comoving number density is kept con- 
stant. The bias evolution in such a model is given by: 



6.2 Small scales 

At smaller scales (r < 1 h~'Mp c) the situation is less complicated 
by the flat large-scale slope. Here lSawangwit et al.l found that a viri- 
alised model gave a better fit to the slightly faster evolution needed 
to fit the small-scale correlation function amplitudes than a comov- 
ing model. But in the present case, the scaling between the AAQ 
and Stripe 82 LRGs in Fig. I13b.b. shows that here the comoving 
model is preferred at small scales over the faster virialised evolu- 
tion. This fits with the more general picture of the Stripe 82 LRGs 
presenting a higher amplitude than expected all the way down to 
the smallest scales. Unfortunately the remaining uncertainty in the 
Stripe 82 LRG luminosity class is still too large to make definitive 
conclusions on this evolution possible. 

6.2.1 HOD Evolution 



where D(z) is the linear growth factor. 



(37) 



Given the uncertainty in ^20 caused by the flat w(d) slope on inter- 
mediate - large scales, we will extend further the studies at small- 
scales, using the HOD model to interpret the small-scale cluster- 
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Figure 17. The predicted SDSS LRG w((9)'s at zimer = 0.35 for the case 
of passively (/no-merge = 1 ) evolving the best-fit HOD of Stripe 82 LRGs 
sample from Zearlier = 1 and the case where central galaxies merging is 
allowed from ZearUer = 1 (/no-raeigc = 0.21), in green dot-dashed line and 
blue long-dashed fine, respectively. The bottom panel shows the ratios of 
the evolved vy(f)'s to the SDSS best-fit, the shaded regions signify the Icr 
uncertainties. 



ing sig nal of the LRGs. Based on the HOD fit at z a 1 , we again 
follow ISawangwit et al H), (and references therein) and test 
long-lived and merging models by comp aring the predict ions of 
these models to the SDSS HOD fi t from lSawangwit et all These 
authors and also IWake et al ] l l2008h found that long-lived models 
were more strongly rejected at small scales (r < 1 h^'Mpc) than at 

intermediate-large scales. 

Again we follow the approach of IWake et al and 
ISawangwit et alj ( 12011) who assumed a fo rm for the conditional 
halo mass function Sheth & TormenI ( l2002h and a sub-Poisson dis- 
tribution for the number of central galaxies in low-redshift haloes 
of mass M such that 



{N,{M)) = 1 . 



C(M) 



where N^.^^ = int(M/M„„„), 



C{M): 







dmN(m,M) <A'c(m)> 



(38) 



(39) 



and N(in, M) is the expression of ISheth & TormenI (|2002|) for the 
condit ional halo mass which generalize those of iLacev & Cold 
( Il993l) . The mean number of satellite galaxies in the low-redshift 
haloes is then given by 



{NAM)} (N,iM)} = S(M) + /„o_„,e,ge [C(M) - <A'c(M)>] , 
where 



5(M): 



Jo 



dm N(m, M) {N^m)) {N,(m)} . 



(40) 



(41) 



and the main parameter is f„o-magi: which is the fraction of un- 
merged low-z satellite galaxies which were high-z cent ral galaxies. 

This model is called 'central-central mergers' in IWake et aTl 

( I2OO8I) . More massive high-z central galaxies are more likely to 
merge with one another or the new central galaxy rather than 
satellite-satellite mergers. 

Setting yi,o-raeiee = 1 mcans that there is no merging of initial 
central galaxies in subsequently merged haloes, so it is similar to 



the passive/long-lived model, /no-meise equals to means that all 
the central galaxies in haloes at high redshift merge to form new 
central and/or satellite galaxies in the low redshift haloes. In the 
analysis below, we use the best-fit HOD model values as estimated 
for scales up to 45' (see Table[3](. 

The /no-mcrgc = 1 case is shown as the w{d) passive model in 
Fig. [T7] and is clearly rejected by the data at < 10' (see lower 
panel). Best-fit HOD predictions of the satellite fraction in the case 
of the passively evolved LRGs from Zeariier = 1 to zimer = 0.35 
is Fj„f = 18.6 ± 2.5% whereas [Sawangwit et al.l measured Fj^, = 
18 ± 1% for a brighter selection of LRGs at Zem-iier = 0.68. We see 
that both these results, for the long-lived model, are significantly 
higher compared to the best-fit SDSS HOD, F,„, = 8.1 ± 1.8%. 
The difference in the number of the satellite galaxies is explained 
as the predicted clustering amplitude at small scales (1-halo term) 
for the passive model, is higher compared to the SDSS HOD fit 
as it is clearly shown in Fig. [17] Higher clustering signal at small 
scales indicates the presence of too many satellite galaxies in the 
low-redshift haloes. 

The wiff) merger model is described by f„o-mcnc = 0.21 as 
presented in Fig.[T7]and clearly fits the data well. For this model the 
satellite fraction at z = 0.35 estimated to be F ^, = 7.29 ± 4.5% and 
is in a good agreement with lSawangwit et al.l Moreover, the best-fit 
HOD model values for the evolved Zeariier = 1 LRGs to zimer = 0.35 
for bias and galaxy number density are b = 2.24 ± 0.24 and jig = 
0.67±0.41 10"*h''Mpc"', respectively. Compared to the SDSS best- 
fit model, with h = 2.08 ±0.05 and iig = 1.3 ±0.4 10"-'h''Mpc-\ the 
number of galaxies at z = 0.35 have been decreased by almost 50% 
due to central-central merging. The evolved linear bias and galaxy 
number density a re consistent with the z = 0.35 best-fit HOD of 
ISawangwit et all at 1 - 1.5 cr level. 

Note that the agreement at large scales in Fig.[T7]is somewhat 
artificial given the underestimation of w(9} by the HOD model in 
Fig. 115b which remains unexplained in the HOD formalism. But 
at these smaller scales the result that the merging model fits bet- 
ter than the long-lived or indeed the virialised clustering model of 
Fig. 113b may be more robust, given the reasonable fit of the HOD 
model at small scales (d < 3') in Fig. 1 15b. 



7 TESTS FOR SYSTEMATIC ERRORS 

In this section we will present an extended series of checks for 
systematic errors that might have affected our clustering analysis, 
with the major issue being the flatter slope at large scales as 
estimated in 35.21 35.31 and 35.41 Tests for possible systematics 
that will be discussed here are: 



• data gradient artefacts, 

• w{6) estimators bias, 

• survey completeness, 

• observational parameters ; such as star density, galactic extinc- 
tion, seeing etc. 



7.1 Data gradients and w(d) estimator bias 

A false clustering signal at large scales can arise from artificial gra- 
dients in the data, as the correlation function is very sensitive to 
such factors. In attempting to explain the behaviour of the observed 
w(d) at large scales, first we divide the LRG sample area in 6 equal 
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Figure 18. m'(9)'s from Landy Szalay, Hamilton and standard estimator of 
the 700 deg"^ LRG sample. For comparison, the averaged »'(6)'s from the 6 
subfields (see text for more detail), are overplotted as measured from each 
estimator. Landy & Szalay and Hamilton estimators, give similar results for 
the average subfields and the full sample measurements, respectively. The 
standard estimator is more biased, at larger scales. 



Figure 19. w{6) results of the 6 equal size subfields (15 X 2.5deg each) 
across Stripe 82, the total area as estimated by using the Landy Szalay esti- 
mator and the averaged clustering signal from the 6 subfields. In the bottom 
panel are displayed the ratios of the w{6) of each subfield compared to the 
total area. 



K 



LRGs 700 deg- 



17.0-17.2 
17.2-17.4 
17.4-17.6 
17.6-17.8 
17.8-18.0 



4894 
11096 
22490 
38659 
53680 



Table 4. K-limited sub-samples used for auto-correlations in Fig. 20. 



subfields in RA. Then the angular correlation function of each sub- 
field has been calculated using the Landy & Szalay, Hamilton and 
the Peebles estimator - the standard estimator. Furthermore, we av- 
erage the w{6) results of the 6 subfields as measured by each esti- 
mator and we compare them with 700 deg"^ LRG w(Q) full sample 
results (see Fig.llSt. 

From these comparisons, it is clear that when we use the 
Landy & Szalay and Hamilton estimators, we do not find any sig- 
nificant difference in the amplitude of the measured w(0) between 
the averaged subfields' or between the full samples' measurements. 
When the averaged w{6) measurements are compared with those 
from the full sample, only a very slightly smaller clustering signal 
in the averaged w(f?)'s is seen, barely visible in Fig. [18] Further- 
more, this is only the amount expected from the integral constraint 
(see ^3]2} on w(6), if the above Landy & Szalay estimate is as- 
sumed to apply in a single sub-field area. The standard estimator is 
known to be subject to larger statistical errors at large scales and 
here the signal is actually stronger when compared with the other 
two estimators. 

Moreover, in Fig. [19] we display the results of the w{Q) mea- 
surements from the 6 subfields individually against the full sample 
measurements as estimated with the Landy & Szalay estimator in 
all cases. Even now we cannot see any major trend through the sub- 
fields' correlation function measurements, except possibly for the 
15" < RA < 30" subfield which has a steeper slope at larger scales. 



700deg LRGs sample o 
17.0 < K < 17.2 □ 
17.2 < K < 17.4 X 
17.4 < K < 17.6 ^ 
17.6 < K < 17.8 □ 
178 < K < 18.0 A 



1.0 10.0 
e/arcmin 



Figure 20. Auto-correlation functions from Landy-Szalay estimator for the 
700 deg^^ LRG K-limited sub-samples from Table 4. Total sample is over- 
plotted for comparison. 



7.2 Magnitude incompleteness 

Another issue that we want to address is how the clustering sig- 
nal can be affected by magnitude incompleteness. The izK colour- 
selection used for the LRGs, applied up to the faintest limits of the 
SDSS-UKIDSS LAS surveys (see ^211 . To account for this, first 
we divide the 700 deg"' LRG sample in 5 A" magnitude bins in the 
range 17 < A" < 18. The number of LRGs in each magnitude bin is 
shown in Table [4] 

Measurements of the angular correlation function from each 
K-hin are shown in Fig. [20] where measurement uncertainties are 
not shown as we are mostly interested in the shape of the w{d) in 
the linear regime. The clustering signal from the A"-magnitude bins 
compared to the full sample do not show any significant difference 
at large scales and follow the full sample w(d) shape. At smaller 
scales we see that the clustering from the brighter samples is higher 
than for the fainter samples, as expected. 
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Figure 21. a) The obsei'ved w{8) of Stripe 82 LRGs (blue dashed line). Stripe 82 star catalague of Ivezic et al. (2007) autocon'elation (green hne) for < 21, 
cross-coirelations of the aforementioned LRGs-stars (black dashed-dot line) and the resulted corrected observed autocorrelation function following Ross et al. 
(2011). We see that there is no difference between the observed LRGs and the corrected »'(6)'s, respectively, b) The number density of the star's up to ; = 21 
from Ivezic et al. (2007) catalogue (blue diamonds) and the 700 deg"" LRG sample (red triangles) across the Stiipe 82. There is a strong gradient in the star 
distribution towards one end of the Stripe 82 at 330 < RA < 340deg or -2 < RA < -Ihr in the abscissa notation. But when we excluded this area from the 
star-LRG cross-correlation, there was no change in the large-scale w{9) signal. 
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Figure 22. The projected number density of Stripe 82 LRGs as a function of the potential observable systematics: stellar' density (Him, ), Galactic extinction 
(At) in the K-band, the K-band seeing (seeic), K-band background median sky flux in counts per pixel and the airmass (air). The en'ors are the standard 
deviation of the measurements within each bin. 



The final tests of the magnitude incompleteness check are via 
the use of brighter colours in the zK selection. We therefore se- 
lected on the basis of brighter magnitudes down to z < 21.2 and 
K < 17.2, in various combinations and re-measured the angular 
coiTelation function. Even with these bright cuts, we did not see 
any change in the excess at large scales. 



7.3 Observational parameters 

The final test to identify a potenti al observational sy stematic effect 
follows the approach described bv lRoss et al.l ( l201ll) . referring pri- 
marily to the area effectively masked by stars with magnitudes sim- 
ilar to the galaxies in the field. We cross-correlated th e 700 rfeg~" 
LRG sample with the Stripe 82 star catalogue from llvezic et al] 
( |2007|) . in 4 magnitude bins, / < 19.5, 20, 20.5, 21. From the mea- 
sured autocorrelation function of stars and the cross-con'elation 
function of stars with LRGs we computed the effect of stellar mask- 
ing on the LRG correlation function using their equations (28) and 
(29). We show these results in Fig. 121b. 



The cross coiTelation results show a very small anticon'elation 
between LRGs and stars for the ; = 19.5 and 20.5 bins. A possible 
explanation for this anticorrelation might be related to the fact that 
we see an increase in the star number density between 330 < RA < 
340 rfeg" (see Fig. |2T]?). Next, we ca lculate the expected w(Q), as 
defined in Eq. 29 of lRoss et al.l ( 1201 ll) . In all cases, there was little 
difference in the expected and observed w{8) of the 700 deg"^ LRG 
sample. We conclude that the effect of stellar masking is essentially 
negligible, less than 1 % of the clustering signal at a; 90' 

The re are other source s of possible systematics as well as star 
masking. IRoss et aT U20llh also checked observational parameters 
such as galactic extinction, sky background, seeing and airmass 
using the cross-correlation technique. The Stripe 82 LRG sample 
is K-limited. Hence, we explore if the above observed parameters 
from the UKIDSS LAS K-band could be sources of systematic er- 
rors at large scales. Fig.l22ldisplavs the number density of Stripe 82 
LRGs and how it is rel ated with ea ch potential observational sys- 
tematic (stars are from llvezic et al.l 2007). From Fig. |22]we see a 
sharp decrease in the number of LRGs with high galactic extinction 
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and poor seeing. The airmass fluctuations are also large compared 
to the error bars. The majority of the LRGs lie within the first few 
bins of galactic extinction, seeing and airmass in Fig. |22] but the 
LRGs in the rest of the bins with higher values could introduce 
syst ematics in the cl ustering signal. 

IHo et al] j2012l) present a method to identify which combina- 
tion of the observed parameters could have the biggest effect on the 
clustering measurements. The authors in this work expressed the 
linear relationship between the potential observational systematic 
and its effect on the observed overdensity of galaxies, through the 
e factor. In Fig. |23fa we show the e, parameters f or each of the the 
observational parameters. The IRoss et al.l ( I2OI ll) cross-correlation 
correction technique requires that e be constant, so we use the best- 
fitting constant value of e as calculated with the lowest chi-square 
fits from field-to-field eiTors. We find that the biggest correction 
in the angular coiTelation function is for the combined seeing, air- 
mass and galactic extinction observational parameters (see Eq. 29 
of Ross et al.). Also, a slightly smaller correction has been found 
for stars, sky background and galactic extinction. In Fig. |23b we 
show the original uncorrected w{6) for the Stripe 82 LRGs, the w(6) 
after applying the combined correction for the seeing, airmass and 
galactic extinction. In the same figure, for comparison we plot the 
best-fit ACDM model as displayed in Fig. [141 So far this correction 
in our w{ff) results is the most important. But still as we can see 
from Fig. |23b. at the range 20 - 80', the amplitude of the angular 
correlation function does not show the expected behaviour of the 
standard model. We have checked for the most common sources 
of systematics in the literature. Our data could still be affected by 
hidden artefacts, a case that future studies might be able to identify, 
but for the moment we will take the connected result in Fig. 123b as 
our best estimate. Note that the HOD fits of 15. 41 were only done 
up to < 60' where there is little change in the form of our w{6) 
result. 



0.6 
0.4 
0.2 
" 0.0 
-0.2 



O starp 
O Ak 
□ seeihg 

sky 
A air 



20 



40 



60 

) (arcmin) 



80 



100 



0.0100 



0.0010 




^ Uncorrected 
Seeing + Airmass + Galactic Ext corrected 
ACDM 



) (arcmin) 



8 TEST FOR NON-GAUSSIANITY 

One possible explanation for the flat slope seen at large scales is 
scale-dependent bias, although this is usually discussed more in the 
context of small-scale clustering. However, scale dependent bias at 
large scales has previously been invoked to explain the discrep- 
ancy between the APM w{ff) results and Q.„, = I CDM models 
teower et al.lll993l) : in this case the scale dependence was caused 
by 'cooperative galaxy formation'. 

Another possibility is that the LRG power spectrum may be 
closer to the primordial power spectrum at higher redshifts . But we 
have seen that the Stripe 82 clustering result are not in line with 
the standard ACDM model. These correlation function results are 
better fitted by a model with Cl„, = 0.2 rather than CI,,, = 0.27 (see 
Fig. I14t . useful at least as an illustration of the size of the LRG 
clustering excess. 

The third possibility is that the z ~ I LRG power spectrum 
may be better explained by scale-dependent bias at large scales due 
to primordial non-Gaussianity in the density fluctuations. The pri- 
mor dial non-Gaussianit y of the local type is parameterised by /^"l^' 
(see lBartolo et al .|2004 for a review) and is expected to contribute a 
l/k^ term to the power spectmm and evolves as a 1 -l-z (see Eg. 142b. 
It is therefore best seen at large-scales and high redshifts. Fig. 1 
of IXia et"^ 1 I2OIOI) shows the potential effect of non-Gaussianity 
on the biased clustering of radio sources with a similar redshift to 
the LRGs discussed here. It can be seen that the non-Gaussianity 



Figure 23. a) (upper) Similar to Ross et al. 201 1 we plot e, the linear factor 
between the potential observational systematic and its effect on the observed 
overdensity of galaxies for stars (puiple diamond), galactic extinction (blue 
diamond), seeing (red squares), sky background (green diamond) and air- 
mass (orange triangle). The solid Unes are the best-fitting constant value of 
£ for each systematic, b) (lower) The w{0) measurement of the Stripe 82 
LRGs without any cross-con'elation connection (black star) and w(6) cor- 
rected for seeing, airmass and galactic extinction combined (red diamond). 
The best-tit ACDM model to the uncorrected measurement is plotted (blue 
line). 



causes a strong positive tail to the coixelation function for 9 > a 
few degrees. 

IXia et al.l ilOld) . foflowing IsTake & Walll l [20o3) . found that 
the NVSS surve y ACE showed a strong positive tail suggesting that 
= 62 + 27. IXia et al.l ( l201lh also inspected the angular corre- 
lation function of the DR6 I x 10* QSO sample and found similar 
results to the radio sources with again an extended correlation func- 
tion being seen implying similar values of /nl (hereafter we shall 
use just /nl to denote ./jjjjl;"') as for the radio sources. This led to 
only slightly weaker constraints than for the radio sources in terms 
of th e value of /nl- 

ISawangwit et alj ( I2OI ih measured the combined angular cor- 
relation function of LRGs at z a 0.35,0.55,0.68 and found that 
although the results were in agreement with ACDM at scales < 
lOOh^'Mpc, at larger scales there was a possible excess, although 
this could still be due to systematics. 

We then proceeded to follow Ixia et al.l and fit /nl models. We 
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Figure 24. a) (upper) The combined congelation function of Sawangwit et 
al (201 1) for the z = 0.35, z = 0.55 and z = 0.68 LRG samples, compared 
to a standard LCDM model (/„; = 0) and models with increasing degrees of 
primordial non-Gaussianity (/„/ = 62, 80). b) (lower) The Stripe-82 ; a 1 
LRG congelation function compared to a standard LCDM model {f„i = 0) 
and models with increasing degrees of primordial non-Gaussianity (/„/ = 
62, 80,100). 



use their relation between the non-Gaussian and Gaussian biases 
(bNG and ba) 



hiaiz) - baiz) - 2(ha(z) - 1),/nl5cc(z)Q'm W ■ 



(42) 



Here 5ec(z) is the critical density for ellipsoidal collapse and 
ayi{k) <xi l/k^ contains the scale and halo mass dependence (see 
IXiaetalJ for more details.) 

We first applied this relation to the case of the NVSS radio 
sources at z a; 0.7. We found that adding the l/k^ term to the stan- 
dard cosmology P(k) caused it to diverge and so we had to apply a 
large-scale cut-off, so that for k < kg then P(k) = 0. This is clearly a 
source of uncertainty in fitting for /nl- Neverth eless, we found th at 
for ^To = 10"'', we could reproduce the results of lXia et all ( 1201 Ol) . 

We then applied the same technique and cut-off to the com- 
bined AAQ. LRG and the Stripe 82 LRG w(e)'s (after applying the 
combined correction for seeing, airmass and galactic extinction as 
estimated in § |7.3t. We first took the value of b^ = 2.08 from the 
halo model fits of ISawangwit et al. l|) and fitted for fi^t. The 



Figure 25. The minimum^- is 5.5 over 1 1 d.o.f and the best-fit parameters 
are /nl = 90 ± 30 (Icr) and M„,„ = 1.26 ± 0.22 x 10"h-'Mo. The best- 
fit M„„„ here is lower than the full HOD fit assuming fyi^ = at 2.2 X 
lO'^'^h-'Mo. 



result is shown in Fig.l24fa. We find that for AAQ LRGs, the results 
for /nl are reasonably compatible with those from the NVSS cata- 
logue with values of /nl = 60- 80 giving a better fit to the data in 
the range 1.5 < d < 6.5deg. 

The prediction from non-Gaussianity is that the large scale 
slope will further flatten with redshift. We therefore compared the 
Stripe 82 LRGs to models with the same /ml values and find no 
inconsistency (see Fig. 124b). Clearly the errors at the largest scales 
are more significant for the Stripe 82 data than for the AAQ. LRGs 
or the NVSS radio sources. However, the predicted flattening of the 
Stripe 82 correlation function at (? a Ideg makes the non-Gaussian 
models more consistent with the data in this smaller angular range 
than the /nl = model. At larger scales the errors are larger and 
the data is therefore more in agreement with the standard model. 

Fig. |25|shows the effect of jointly fitting /^vl on the minimum 
halo mass, M,„„,, in the HOD model. The best fit model now gives 
M,„,„ = (1.26 ± 0.22) X lO'^/i-'Mo and ffjL = 90 ± 30, lower than 
then the M„„„ = 2.2 x IO'^A^'Mq value when /,/ = is assumed in 
the full HOD fit. 

We should say that rather than detections of non-Gaussianity, 
the present AACl and Stripe 82 LRG results should be more re- 
garded as upper limits on non-Gaussianity. Large-scale angular cor- 
relation function results are still susceptible to large-scale gradi- 
ents and even though there is no direct evidence for these in the 
AAQ or Stripe 82 samples, there is still the possibility that these 
exist in the data. On the other hand, the classic test for the real- 
ity of a correlation function feature is that it scales correctly with 
depth and at least the SDSS and Stripe 82 LRG correlation func- 
tions in Figs. I24fa.b lo ok like they do so. It w ill be interesting to 



see if as QSO surveys I Sawangwit et al .120 12h and z ~ 3 LBG sur- 
veys teielbv et al.l20ia grow, whether the correlation functions at 
higher redshift also show an increased slope flattening as predicted 
for the non-Gaussian models. 

The other uncertainty that has arisen is in the non-Gaussian 
model itself where we have found that there is a rather strong de- 
pendence on a small-scale cut-off, kg. Other authors have made 
some reference to this problem but only implicitly. It will be in- 
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teresting to see if more accurate models for non-Gaussianity can 
numerically predict this cut-off from first principles. 



9 SUMMARY AND CONCLUSIONS 

We have measured w{0) for a; 130000 colour selected galaxies in 
Stripe 82 exploiting SDSS DR7 ; + z bands and UKIDSS LAS 
K photometry. We used the cross-correlation technique of New- 
man (2008) to establish that the average redshift of the LRGs is 
z ~ \. This sample therefore probes higher redshifts than the pre- 
vious SDSS LRG samples of ^awangwit et al. (2011^). We have 
established that a sample with sky density a TOOdeg"^ has a 
comparable space dens ity to the z ~ 0.68 AAQ. LRG sample of 
ISawangwit et alj ( l201ll) . However, this is only an approximate cor- 
respondence which makes evolutionary comparisons between the 
redshifts more tricky. What is clear is that the z~ I LRGs generally 
have a relatively high clustering amplitude. Compared to the AAQ. 
LRG w(e) scaled to the depth of the Stripe 82 LRGs, the Stripe 82 
wifl) is higher at all scales, even those below < 1 h^'Mpc. Thus 
at intermediate scales, the z « 1 LRGs are not only more clus- 
tered than predicted by the long-lived evolutionary model, they are 
also more clustered than the comoving model. At small separations 
(< 1 h^'Mpc) the correlation function amplitude is again somewhat 
higher than the AAQ. results scaled by the previously preferred sta- 
ble clustering model. The Stripe 82 w(Q) also shows a very flat slope 
at large scales which means that the ACDM linear model has be- 
come a poorer fit than at lower redshift. 

Partly to look for an explanation for the flat large-scale slope, 
we then fitted a HOD model to the Stripe 82 w(e). The best fit 
parameters were M,„,„ = 2.19 ± 0.63 x 10" h^'Mo, Mi = 21.9 ± 
5.6x 10" h-'Mo, fci,„ = 2.81 ±0.18, M,g = 3.3±0.6x 10" h-'Me, 
Fsa, = 3.17 ± 0.08% and ;ig = 0.8 ± 0.3 x lO--* h^Mpc^^ The high 
amplitude of the correlation function clearly pushes the halo masses 
up and the space densities down. The lowest chi-square fits were 
found when large scales were excluded but the reduced chi-squares 
were still in the range 2.3-3.6. This is actually an improvement over 
the lower redshift samples but this is certainly due to the larger 
errors on the Stripe 82 data. We conclude that it is not possible to 
find an explanation for the flat slope in the Stripe 82 w{6) on the 
basis of the HOD model. 

We also then studied the evolution of the HOD between z = 1 
and z = 0.35. Similar to ISawangwit et al.l ( l201lb . we concluded 
that a pure passive model with a low merger rate might produce too 
steep a w{6) slope at small scales (< 1 h^'Mpc). In this case, we 
have already noted that the small scale amplitude may also be too 
high for a passive model with stable clustering. 

We have looked for an explanation of the flat slope in terms 
of systematics by cross correlating the Stripe 82 LRG sample with 
stellar density, airmass, seei ng, sky background and galactic extinc- 
tion and used the method of IRoss et al.l j201lh to correct our w(d). 
Even the combined correction for seeing, airmass and galactic ex- 
tinction only produced a small change in w{6) at large scales. 

We conclude that the high amplitude and flat slope of the 
Stripe 82 LRGs w(ff) may have a significant contributions from the 
uncertainty in the comparison between AAQ. and Stripe 82 LRG lu- 
minosities. However, this leaves a similar contribution from a new 
and unknown source. We have discussed large-scale, primordial, 
non-Gaussianity as one possibility for the source of this large-scale 
excess. We have suggested that the evidence from the AAQ sam- 
ple itself for an excess at even larger scales may fit in with the 
behaviour expected from non-Gaussianity over this redshift range. 



In this case we returned to the fitting of halo masses including the 
non-Gaussian component and found that the best fit M,,,-,,, decreased 
from 2.2 x 10"Mo to 1.3 x 10"Mq. More importantly, if the Stripe 
82 large-scale w{6) excess proves reliable and not due to systemat- 
ics, then we have made a significant detection of non-Gaussianity in 
the z a 1 LRG distribution with an estimated local non-Gaussianity 
parameter estimate of /,1,'J^''' = 90 ± 30. This represents a 3cr detec- 
tion at a level comparable to the presen t upper limit fro m WMAP 
CMB measurements of f^^^ = 32 ± 21 jKomatsj|201Cl) . 
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